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Abstract. We derive an extension of the standard time dependent WKB theory 
which can be applied to propagate coherent states and other strongly localised 
states for long times. It allows in particular to give a uniform description of the 
transformation from a localised coherent state to a delocalised Lagrangian state 
which takes place at the Ehrenfest time. The main new ingredient is a metaplectic 
operator which is used to modify the initial state in a way that standard time 
dependent WKB can then be applied for the propagation. 

We give a detailed analysis of the phase space geometry underlying this 
construction and use this to determine the range of validity of the new method. 
Several examples are used to illustrate and test the scheme and two applications 
are discussed: (i) For scattering of a wave packet on a barrier near the critical 
energy we can derive uniform approximations for the transition from reflection 
to transmission, (ii) A wave packet propagated along a hyperbolic trajectory 
becomes a Lagrangian state associated with the unstable manifold at the Ehrenfest 
time, this is illustrated with the kicked harmonic oscillator. 
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1. Introduction 

The semiclassical propagation of wave packets is of fundamental importance in many 
applications of quantum mechanics and has therefore been studied extensively in the 
literature [l}|9] . Time dependent WKB approximations apply to extended initial states 
of the form 

V'(x) = A(2;)e«^(^) (1) 

where the phase function S{x) is real valued and smooth and the amplitude A{x) 
is smooth and non-oscillating. The time evolution of such states can be expressed 
in terms of transport along a family of classical trajectories determined by initial 
positions x and initial momenta p — VS{x), see, e.g., [2|[5]. 

A different class of initial states is given by coherent states, see, e.g., [3][7j, which 
are strongly localised and therefore their propagation can be described up to certain 
times using only one trajectory and the linearised motion around it. A Gaussian 
coherent state is a function of the form 

s ^ (det Im BY^^ . [p.(^_^)+(^_^).B(^_^)/2] 
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where Z — (p, q) £ R^" is a set of parameters which determine where the state is 
locahsed in phase space and B is a, complex symmetric n x n matrix with ImB > 0. 
The latter condition ensures that the state is localised around x = q. The roles of B 
and Z become clearer if we look at the Wigner function of such a state, which is a 
Gaussian of the form 

W{z) = (i);.e-(^-)-(^-)/'' , (3) 
with phase space coordinates z G E," x R" and where G is a positive definite symmetric 



and symplectic matrix determined by B. As was described by Hepp and Heller 10 - 13 
the time evolution of a coherent state can in the semiclassical limit be described by 
using just the classical trajectory Z{t) = {p{t),q{t)) through Z and the linearised 
flow around it which is a time dependent 2n x 2n symplectic matrix S{t). This simple 
description of the semiclassical propagation of coherent states made them a very useful 
tool in many applications, in particular in chemistry due to the work of Heller and 
coworkers and their appearance in initial value representations like the Herman-Kluk 
propagator [3||4)[l3}|l8] 

We will be using time dependent WKB methods developed for states of the form 
([T]) to understand the propagation of the states ([2| for long times. In order to identify 
the relevant time scales it is useful to consider the Wigner function of the time evolved 
coherent state which in leading order is given by 

lYU ^) ^ 1 ^-{z-Z(t))-G{t)(z-Z(t))/h 

where G{t) = {t)GS (t) and S{t) is the linearised flow around the central trajectory 
Z{t). The basic idea leading to this approximation is that since a coherent state 
is localised around a point in phase space one can approximate the Hamiltonian by 
its quadratic Taylor expansion around the classical trajectory Z{t). The resulting 
Schrodinger equation can then be solved explicitly. This approximation can only be 
expected to be accurate as long as the propagated state stays localised, and from 
^ we see that this is the case as long as \min[G{t)]/h ^ 1, where we denote by 
^min/max[G{t)] the Smallest or largest eigenvalue of the matrix G{t) = 5*(t)G5(t), 
respectively. Since G{t) is symmetric and positive we have Xmin[G{t)] > 0, and 
that G{t) is symplectic implies that l/A.,„.m[G(i)] = Xmax[G{t)], so the localisation 
condition becomes 

XmaAG{t)]h ^ 1 . (5) 

In order to relate this more explicitly to the properties of the classical dynamics we 
use the estimate Xmax[G{t)] = ||G(i)|| < G||5(i)|p and so the condition becomes 
\\S{t)\\Vh < 1. Here we use for a matrix A the norm induced by the Euclidean 
norm ||a:|| = ^yx ■ x, i.e., \\A\\ := sup||j.||^]^ IMF' time at which the width of a 

propagated coherent state reaches order one is called the Ehrenfest time TE{fi), and 
by the previous discussion we see that it is determined by the condition 

\\S{TE{h))\\Vh^i . (6) 

Hence the Ehrenfest time depends on the dynamical properties of the classical system 
near the trajectory Z{t). If the trajectory Z(t) is hyperbolic, with largest Liapunov 
exponent A > 0, then ||S'(i)|| ^ e"^* and so 
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on the other hand side, if the dynamics is integrable we typically have ~ t for 

large t and then 

TEih) = . (8) 

Ehrenfest showed in his classical paper [l9] how one recovers classical mechanics 
from quantum dynamics by considering expectation values in propagated localised 
states. His construction breaks down if the state becomes delocalised, hence the name 
Ehrenfest time for the time scale at which this happens. In more recent years the 
behaviour of time evolved coherent states for long times has attracted the attention of 
mathematicians and rigorous estimates for long times have been derived, in particular 
in 20 21 it was shown that the remainder terms stays small for times satisfying 
|t| < ^Te in the hyperbolic case. One should emphasise however that the breakdown 
of Ehrenfest 's argument at the Ehrenfest time does not mean that the quantum to 
classical correspondence breaks down, it only means that the approximate propagation 
of a wave packet based on one central trajectory breaks down. If one includes more 
trajectories one can use semiclassical propagation well beyond the Ehrenfest time as 
has been demonstrated in p2{[23] , although from a mathematical point of view this is 
still a challenging problem. 

The principal aim of this work is to present a propagation scheme for coherent 
states which works at times of the order of the Ehrenfest time t ^ Te, and is able to 
describe the transition from a localised state ^ to an extended state ([T]) at the 
Ehrenfest time in a uniform way. In a previous paper [24| it was demonstrated 
that at the Ehrenfest time a coherent state which is transported along a hyperbolic 
trajectory becomes effectively a WKB state of the type ([T]) associated with the unstable 
manifold of that trajectory. Here we will give the theoretical foundation for that claim 
presenting a new propagation scheme which is a combination of a time dependent 
WKB approximation and a metaplectic correction. 

The qualitative change in the nature of the state at the Ehrenfest time was 



described in |25 for the cat map and in 26 for the propagation of coherent states 



centred on unstable fixed points in one-dimensional multiple well potentials. In 
contrast to the present approach that work was not based on direct propagation of 
the state in the position representation but on the Wigner function and the Egorov 



theorem from 27 . That means in particular that we can give a much more detailed 
description of the state at and beyond the Ehrenfest time. 

The plan of the paper is as follows. In Section [2] we recall time dependent 
WKB theory for the propagation of states of the form (iTl and rewrite it in an exact 
form where we explicitly separate the classical transport and the quantum dispersion. 
The main idea of this paper is then developed in Section |3] where we show how to 
approximate the action of the dispersive part on a coherent state using a simple 
metaplectic operator and combine this with time dependent WKB approximation 
to obtain a method which allows to describe the propagation of coherent states at 
and beyond the Ehrenfest time. In Section [4] we develop the phase space geometry 
underlying the metaplectically extended WKB method. This allows us to determine 
its range of validity, in particular its dependence on dynamical properties of the 
classical system. We then develop the geometric picture even further in Section [s] 
to allow more general Hamiltonians and the inclusion of caustics after the Ehrenfest 
time, the price we have to pay is that remainder estimates become less explicit 
now. The last two sections are devoted to examples. In Section [6] we consider a 
couple of simple explicitly or almost explicitly solvable examples which demonstrate 
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in some detail how our method works in practice. We show in particular that for 
a potential barrier we can describe the transition from transmission to reflection 
near the critical energy in a uniform way. In Section [7] we reconsider the situation 
from 24 and demonstrate how our method reproduce the results. As a particular 
application this shows using 28 that in a strongly chaotic system coherent states 
become semiclassically equidistributed beyond the Ehrenfest time. In the end we 
summarise our results in the conclusions. In the appendix we collect some more 
technical results from semiclassical analysis we use in the main body of the work. 



2. Time dependent WKB 

Our method is an extension of the well known time dependent WKB method for 
real valued phase functions. So we start by recalling this method, and to keep the 
discussion simple we restrict ourselves first to the case of the standard Schrodinger 
equation with potential V in R", which in appropriate units reads 



ihdti/j 



-At/j + FV' 



(9) 



in Section [5] we will discuss a more general setting. We want to solve this equation for 
an initial state of the form tp = AqBTi^", where Sq is real valued. Inserting the usual 
WKB ansatz 

= A(^,a;)e^^(*■^^ (10) 

with S{t,x) real valued, into the Schrodinger equation gives an expression which we 
separate, following the standard treatment into the following two equations 



dtS+-\VS\ 



V = 



idtA + iVS ■ VA + h{AS)A = ^AA . 



(11) 
(12) 



If A does not depend on h, then this is a splitting into different powers of h. 



Equation (111 is the Hamilton- Jacobi equation and its solutions are described using 



the propagation of the Lagrangian manifolds defined by p VS'(i, q), i.e. 

At := {{WS{t,x),x)} c R" X K," , 



(13) 



which are transported by the Hamiltonian flow associated with the classical 
Hamiltonian g) = ^p'^+V{q). We will analyze this further in Section|4| The phase 
function S{t^x) often exists only on a finite time interval and develops singularities if 
caustics appear, then one has to use a refined Ansatz. Because of the association with 



Lagrangian submanifolds WKB states of the form ( 10 1 are often called Lagrangian 
states in the mathematical literature 



29 30 



If we ignore in ( 12 ) the term of order h on the right hand side, then this is a 



transport equation which describes how the amplitude A is transported along the 
Lagrangian manifold Aj. Let us define the transport operator T{t) as the solution to 
the transport equation 



idtT{t) = 



iVS"- V 



a'-as 



T{t) , T(0) = / 



(14) 
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where the operator in brackets on the right hand side is self-adjoint, hence T{t) is 
unitary. Then if we make for A{t) an ansatz 

A{t) = T{t)D{t)Ao , (15) 



with an operator D{t) which is to be determined, and insert this into (12), then we 
obtain for D{t) the equation 

id,Dit) = -^Ait)Dit) , (16) 

where 

A{t):^T*{t)AT{t) , (17) 
and with initial condition D{0) — I. Since A{t) is self-adjoint, D{t) is unitary as well. 



From a more detailed analysis of the operator T{t) in Section 4.2 we will learn that 
A{t) is a generalised Laplacian of the form 

A{t) ^^a,j{t,x)d,dj , (18) 



see ( 49 1 . Collecting all terms we have a solution to the original Schrodinger equation 
of the form 

m = [T{t)D{t)Ao]e^-^^''> . (19) 

So the time evolution is described by three parts, (i) propagation of the Lagrangian 
manifold At, (ii) classical transport of the amplitude by the unitary operator T{t) and 
(iii) quantum dispersion described by D{t). Note that both T{t) and D{t) depend on 
the initial state via the initial phase function 5*0. 



The main condition we need in order that ( 19 ) holds on a time interval [0, T] is 



that the projection of Ai to position space is smooth for all t e [OjT], i.e., that At 



does not develop caustics. Under this condition the representation ( 19 1 is exact, and 
S(t) and T(t) are both determined by transport along classical trajectories in phase 
space, the only contribution not yet linked to classical dynamics is the dispersive part 
D{t). We get the usual WKB result by neglecting the dispersion, i.e., by replacing 
D(t) by the identity /. To see what kind of error we make by doing this we can use 
Duhamel's principle which we will state in general form for later use as well [20j. Let 
H{t) and Hi{t) be two time dependent self-adjoint operators, and U{t) and Ui{t) the 
time evolution operators generated by them with initial conditions C/(0) = t^i(O) — I, 
then ^ 

U{t)-Ui{t) = -\ I U{t)U*{s)[H{s)~ Hi{s)]Ui{s)d8 . (20) 



This formula allows to estimate how close the time evolutions generated by two 

2 



different Hamiltonians are to each other. If we choose H{t) = -^A(t) and i? i = 



then U{t) = D{t) and Ui{t) = /, and (20) gives 



D{t)Ao =Ao + - D{t)D*{t')A{t')Ao dt' . (21) 
^ Jo 

Since D is unitary we obtain then directly 

\\Dit)Ao - AoW < ^ j\\A{t')Ao\\ dt' . (22) 
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So if the second order derivatives of Aq are bounded and the coefficients of A{t) are 
not growing with t, then the solution to Schrodinger's equation satisfies. 

V^W = (r(<)Ao)eH^(*)+0(|i|;i) (23) 

which is the standard time dependent WKB result. Furthermore Duhamel's formula 
can be iterated and gives an expansion with error term of order 0{{h\t\)'^) if the the 
first 2N derivatives of are bounded and the coefficients of A{t) don't grow with 
t, (20 . 



3. Metaplectic extension of WKB 

We would like to apply (19) to a coherent state, i.e., a state of the form ([2|. With 
Soix) ^p-ix-q) and Ao{x) = (7r;i)-"/4g-5L(:r-g).B(x-9) ^j^jg ^^^^^ ^^le form we 



can use in (19), but now 

\\A{t)Ao\\ ^ l/h , (24) 



so (22) does not allow us to use the standard WKB approximation D{t) « /. 

The way to solve this problem is to approximate the action of D{t) on Aq not by 
the identity, as in the WKB method, but to borrow from the standard propagation of 
coherent states and approximate the generator of D{t) by its Taylor expansion around 
the centre of Aq up to second order. Since D{t) acts on a state which is concentrated 
in position space at x — q and in momentum space at p = this means that we freeze 
the coefficients of A(t) at x = q and approximate D{t) by the operator generated by 
it. 

To formalise this idea we introduce for q G H" the unitary operator Lq by 

iLqa){x) h-''/^a{{x - q)/Vh) , (25) 

which shifts by q and then rescales in h, so that {Lqa){x) is concentrated around x ^ q, 
e.g., if a{x) = 7r~"/4e~^'^^/^ then — LqO is the amplitude of the coherent state ^ 



considered above. Now inserting A ~ LqO into ( 16 ) gives for a the equation 

idta = hLqA{t)Lqa , (26) 
and we will approximate the operator on the right hand side by 

Aq{t) := \unhL*qA{t)Lq . (27) 



By (18) one finds that Aq{t) — '^aij{t,q)didj and 

Qq{t) := hL*qA{t)Lq - Aq{t) = Ot{Vh) (28) 

in the sense that Qq{t)a{x) = Ot{V^) if a is smooth and has bounded derivatives. So 
the operator Aq{t) is indeed obtained from A(t) by freezing the coefficients at x — q 



and furthermore discarding the first order terms. In (52) we will obtain more explicit 
bounds on ( p8| in terms of the classical flow. 

Therefore if we define the operator Dq (t) by 

idtDqit) = ~^Aq{t)Dqit) , Z?, (0) = / , (29) 
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then we expect that D{t)Lqa ~ LqDq{t)a holds. To check this we use agam Duhamel's 
principle (l20|) 



D{t)Lqa = LqDq{t)a+- D{t)D*{s)LqQq{s)Dq{s)ads (30) 
^ Jo 



and, since D(t) and Lq are unitary, we get from (281 

\\D{t)Lqa-LqDq{t)a\\< f \\Qq{s)Dqis)a\\ds = Ot{Vh) 

Jo 



(31) 



The t-dependence in the remainder term will be governed by the behaviour of Qq(t). 
In the next section we give conditions on the initial state and on the classical dynamics 



under which the coefficients of Qq{t) stay bounded, see (51) and (52) 

To summarise our results so far, if ipQ = AoeTi^° where 
around q, then the time evolved state is given by 

^(t) = [T{t)LqDq{t)a]ei^^'^ +OtiVh) , 



LqO is concentrated 



(32) 



so we can use standard time dependent WKB to propagate coherent states if we 
include the additional operator Dq(t). This operator has a quadratic generator and 
is therefore a metaplectic operator, see [3 31 32 , and so we call it the metaplectic 
correction to the standard time dependent WKB method. In order to understand 
the behaviour of Dq{t) in some more detail we have to look at the propagation of 
the Lagrangian manifold A by the classical dynamics, and the associated geometric 
interpretation of our extended time dependent WKB scheme. This will be the subject 
of the next section. 

We finally note that it is sometimes useful to work with Mq{t) := LqDq{t)L*, 
because then LqDq{t)a — Mq{t)Lqa = Mq{t)Ao and we can allow for Aq to be of a 
more general form than the one induced by the specific scaling with Lq, as long as it 
is concentrated around x = q. Notice that Mq satisfies the equation 



ihdtMqit) = --Aqit)Mq{t) , M,(0) 



(33) 



and the time evolution of an initial state ipo = Aoen^°, where Ao{x) is concentrated 
around x = q, can be approximated by 



^P{t) = [T{t)Mq{t)Ao]cr.sW + Ot{Vh) 



(34) 



4. Geometric interpretation 

In order to understand the scope and the accuracy of the metaplectic extension of time 
dependent WKB theory we outlined in the previous section, we have to understand 
some of the geometry underlying it. We will see that in particular the time dependence 
of the remainder term is governed by the way in which the classical flow transports the 
initial Lagrangian manifold A. Furthermore, a proper understanding of the geometry 
will allow to extend the scheme beyond caustics. 
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4.I. Transport of the Lagrangian manifold 



Let us first return to the solutions of the Hamilton Jacobi equation (11), which we 
will discuss for a general Hamiltonian -ff a;), 

dtS{t,x) + H{\IS{t,x),x)=Q . (35) 

As is well known, there are two main ingredients involved in the solution, see e.g.. 



|5[|29 30 . The first is the transport of the initial Lagrangian manifold 

A := {{VS{x),x) : a; e [/} , (36) 

where U C R" is an open set which contains the support of the amplitude A{x). 
Let us denote by <i>'(^,a;), where (C,x) e R" x R", the Hamiltonian flow generated 
by if, i.e., <I>*(^, x) = (p{t),q{t)) where {p{t),q{t)) are the solutions to Hamiltons 
equation p = —VqH{p, q) and q = \7pH{p, q) with initial conditions p{t = 0) = ^ and 
q{t = 0) = X, respectively. Using this flow we can transport the Lagrangian manifold 
A and get a family of Lagrangian manifolds 

At = $*(A) = {$*(VS'(a;),a;) : a; e C/} . (37) 

We will call the initial Lagrangian manifold A non- contracting (with respect to 
the flow $*) if there exists a C > such that for alH > 

||d$*|A(z)|| > C , for ah zeA. (38) 

Here d$*|A(2) : T^A — ^ r$t(^-)At is the restriction of the linearised flow a,t z — {p,q) 
to A which is given in local coordinates by the matrix of the derivatives of the 
components of $* with respect to the coordinates on A. This condition means that, 
roughly speaking, trajectories starting nearby on A do not coalesce into each other. 
An example for a non-contracting submanifold is the unstable manifold of a hyperbolic 
trajectory, and more generally, if a system is hyperbolic then any submanifold which 
is transversal to the stable foliation is non-contracting. Stable manifolds are then of 
course examples of manifolds which are contracting, and so not non-contracting. In 
an integrable system any manifold which lies in the regular part of the foliation of 
phase space into invariant tori is non-contracting. 

The second ingredient needed to determine S(t, x) is the projection of A^ to 
position space along the momentum directions, i.e., we take the projection tt : 
Rp X R^' —?' Wg deflned by tt{p, q) = q and restrict it to At. If this map, 

^A, : At ^ R" , (39) 

has no singularities, i.e., 7r(At) = Ut does not contain any caustics of At, then there 
exist, at least locally, a phase function S(t, x) such that 

At = {{VS{t,q).,q) : qeUt} , (40) 

for some open set Ut G R". By assumption ttaq : Aq — > [/q C R" is smooth so we can 
invert it, with inverse given by 7r^^^(a;) = {'VS{0,x),x), and hence we can form 



0a (t, a;) := tta, ^' {'Kf.^ix)) 



(41) 
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which is a map from U ^ Ut, and if vtaj is smooth, it is a smooth and invertible 
map. But notice that (f>\{t) is not a flow. Analogously to the above notion of non- 



contracting, (38 1, we say that 0a(^) is non- contracting on U d R" if there exist a 



C > such that 

W/^{t,x)\\>C , foraU xeU, (42) 

independent of t > 0, where 4>'j^{t,x) denotes the matrix of first order derivatives of 
the components of (l)\(t,x). A necessary condition for (jj^ to be non-contracting is 
that Aq is non-contracting, but in addition we need that the projection tt : At — >■ R" 
is not singular. This implies that we have to stay away from caustics. 

4-2. The transport operator 

Another description of the map 4>A{t, x) is as follows: take the trajectory {p{t), q{t)) — 
$*(^,x) which starts &i t ~ Q bX q{t = Q) — x with initial momentum ^ = VS{x), 
then the map (j)p^{t^x) is the q component of this trajectory, (j)A{t,x) = q(t). Since by 
Hamilton's equations q — \'pH{p,q) and on Aj we have p — ViS'(t, g), we find that 
4)A(t,x) is the unique solution to 

^^^^1^ = VpHiVS{t, 0a(<, x)),Mt, ^)) , with 0a(O, x)^x. (43) 

Notice that we allowed here for a general Hamiltonian, and not only one of the simple 
form kinetic plus potential energy. 



We will now show that the transport operator T{t) defined by (14) can be written 
using this map as 

{Tit)A){x) = [det<l>'j,itAl\t,x))]-'/^A{<j,^\t,x)) . (44) 



By comparing this formula with Appendix A. 2 and using (43 ) we see that the operator 



in (44 1 satisfies 

ihdtTit) = k{t)T(t) , (45) 
where K{t) is the Weyl quantization of the function 

K{t-p,q)=VpH{VS{t,q),q)-p . (46) 



In particular we see that for the special case H = ^p + V{q) we have by Appendix 



1 



WS{t,x)W +^AS{t,x) 



(47) 

and comparing with the transport equation (|12|) proves our claim. 



4-3. The metaplectic correction 

We can use this observation to give a more detailed description of the operator 

--Ait)^^-T*{t)AT{t) . (48) 

This is an application of Egorov's theorem and the technical details are given in 
[Appendix A.3.2[ where we show that the Weyl symbol of this operator is given by 

6H{t, e, x)^l^- Ait, .x)e Tr[(0X')r (i, x)]' , (49) 
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where A{t^ x) is a symmetric matrix given by 



(50) 



and {(j^A )i matrix of second derivatives of the I'th component of (f>~j^ . Bounds on 

the coefficients of the operator A(i) play an important role in estimating the accuracy 
of the time dependent WKB propagation and the metaplectic extension, because 
A{t) appears in the error terms (22) and (31). This is where the non-contraction 
condition we introduced above becomes important, by (42 1 we have that if (j)K{t) is 
non-contracting then 

\\A{t,x)\\<C, (51) 

hence the operator A(t) has bounded coefficients. We can now give as well a more 
explicit description of the operator Aq{t). A short calculation shows that for a phase 
space function (^, x) which is quadratic in the momentum ^ we have HL*HLq = Hq 
with Hq{^,x) = H{^,q + y/hx) and so Hq{(,,x) — Hq{^,q) +0{y/h). Applying this to 
SH{t,^,x) = l^A{t,x)^ + 0{tf) gives that 



A,(t) = V, •^(i,(Z)V, , and hL*„A{t)Lq = Aq{t) + Ot{Vh) 



(52) 



If 0a(O is non-contracting then the remainder is bounded uniformly in time. 

Since we have now an explicit expression for Aq{t) we can compute the action of 
the metaplectic operator Dq{t), defined by (29), on a function a. Set 



Ct 



A{s, q) ds 



and let a(^) = / e ^^^a{x) dx be the Fourier transform of a, then 

1 



{Dq{t)a)ix) 



(27r)" 



e 2 



a(Oe'-«dC 



(53) 



(54) 



Similarly we find for the rescaled operator Alq (t) that integrating ( 33 ) gives 

{Mq{t)A)ix) = J e-H«-c^*«i,(^)ei-?d^ , (55) 

where An{£_) = J e~Jt^^A{x)dx. The operators Dq{t) and Mq{t) act as Fourier 
multipliers with Gaussian functions. For an initial amplitude of the form A{x) = 
(7rfi,)~"/''e~l^l /'■^''^ one finds in particular 



Mq{t)A{x) 



^dei{I + iCt) 



(56) 



The role of the initial manifold 

The localised initial states we consider are of the form -0 = {Lqa)eTi^ , and the 
phase function determines the initial Lagrangian manifold which is crucial for the 
semiclassical propagation scheme we presented. The question we want to consider 
now is if the state V' determines the Lagrangian manifold uniquely, or in other words, 
if there might be other functions a and S such that = {Lqa)eT^^ . In that case we 
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would have some freedom in the choice of the initial Lagrangian manifold we use for 
propagation. 

The condition {Lqa)e'^^ — (L^a)e^'^ gives after multiplication by e^^-^ and 
application of L* that 

h{x) = a(a;)es[S('?+^^)-^(9+^^)l (57) 
and by Taylor expansion we see that 

'-[S{q + Vhx) - ~S{q + Vhx)] = '-[S{q)^~S{q)] + ^[VS{q)-VS{q)]x + m{Kq,x) , 

(58) 

where R{h, q, x) is a smooth real valued function. The first term on the right hand 
side gives just a constant phase factor, so if 

V5(q) = V5(<z) (59) 

then d defines a nice smooth function. But this conditions means that the two 
Lagrangian submanifolds A and A, generated by S and S, respectively, intersect at 
(Pj q) = (VS'(g), g). We therefore conclude that in order to propagate a state localised 
around a phase space point {p, q) we can use any Lagrangian manifold through that 
point which is non-contracting. This freedom of choice can be used to select for 
instance initial manifolds for which the propagation becomes particularly easy, e.g., 
some for which no caustics develop. 

In the case that the trajectory through (p, q) is hyperbolic, any Lagrangian 
submanifold which is transversal to the stable directions is non-contracting, and hence 
can be used for propagation. Although it seems most natural to take the unstable 
manifold, any other manifold which is sufficiently transversal to stable direction will 
converge exponentially fast to the unstable manifold and hence should work with 
comparable efficiency. We test this with the example in Section [7j 



4-5. Relation to standard coherent state propagation and time scales 



For times which are short compared to the Ehrenfest time our approach should 
reproduce the standard results on coherent state propagation. For the general class of 
initial amplitudes of the form (25) the propagation is reviewed in [33,. A propagated 
coherent state is of the form 



'4>{t,x) 



1 

W4 



ah t, 



(60) 



where {p(t),q{t)) is a classical trajectory, l{t) an action type phase which includes as 
well Maslov terms, and ani^t, x) = ao(i, x) + y/hai{t, x) + ■ ■ ■ with leading term given 
as the solution to 



idtaoit.x) 



-I- 



1 



x-V"{q{t))x 



ao(t,x) , ao(0,a;) = oq 



(61) 



This means that the centre of the state is propagated along the trajectory {p{t), q{t)) 
and its shape changes according to the linearised dynamics around the centre. 
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If we expand S{t, x) around x — q{t) up to second order then we find that we can 
write our approximation ( 32 ) in the form ( 60 ) with 

an{t,x) = 6(^,a;)e5-s"(t^«(*))-+o(^/S) , where b{t,x) = [L;i^,^T{t)LMt)a]ix) (62) 
Now using the operator identities 

VhL;^,)^L^^t) = V , (64) 
Ll(^t^VSit,x)L^^t) ^VS{t,q{t)) + VhS"{t,qit))x + 0{h) , (65) 
L;^,)T{t)L,A,L;T*{t)L,^,) = A + 0{Vh) , (66) 

which all follow by direct computation, we obtain with ^S(t,q{t)) = p[t) = q(t) that 



idtKt, x) = -[ix- S"'V + ^ A5 + ^ A j b{t, x) + 0{Vh) 



(67) 



where S and its derivatives are all evaluated at a; = 9(i)- If we furthermore use 
) = [Ab{t, x) + 2iVb{t, x) ■ S"x + (-|S"'a;|2 + iA5)5(i, x)]e5^-^"^ and 



x-S"x\ _ 



A{b{t,x)e2 

combine this with (|67|) we find 
idtah{t,x) 



- -A - -X- S''x- -\S''x\^ 



ah{t,x) + o{Vn) 



(68) 



Finally from the Hamilton-Jacobi equation ( 11 ) we get —x ■ S"x — \S"x\ /2 — x ■ V"x, 
hence the leading term aQ{t,x) satisfies (61). 

This was a formal computation which showed that our result reproduces the 
previously existing results, but it didn't gave much insight where the standard 
approximation might break down. To see this more clearly let us just look at the 
amplitude T{t)Lqb with 6 = DqO. By Taylor expansion around x = qit) = (j>A{t, q) we 
have 4>~^^{x) — g « [<j)'f^{t,q{t))]^^{x — q{t)) and hence 



T{t)Lgb{x) 



1 



b{[VhMt.qm-H^-qm 



(69) 



hi^dct(j)'Jt,q{t)) 

From this expression we see that the state will be localised around x — q{t) if 

VhU'^{t,q{t))\\<^l , (70) 

and since 4>i^{t,x) is a projection of the flow from phase space this is satisfied if 
t <^ Te{K). On the other hand side, if ^/h\\(l)'f^{t,q{t))\\ ~ 1, then the amplitude is no 
longer localised around a point, and the state becomes a Lagrangian, or WKB state. 



5. A generalised propagation scheme and caustics 

So far we have avoided the discussion of caustics, but although there are situations 
where no caustics occur, in many interesting situations we expect caustics to develop. 
A caustic is the image in position space R" of the singularities of the projection 
TTAj : At — >• R", that means at a caustic the tangent plane T^At to At at a 
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point z — {p,q) E At becomes vertical in the sense that T^At D Vz ^ {0} where 
Vz = {(Cj?) ; C G K-"} denotes the vertical subspace. In a neighbourhood of a caustic 
the manifold At can no longer be represented by a generating function S{t, x) as in 
( [40| , and so the simple time dependent WKB propagation we presented above can no 
longer work. This problem is usually solved by switching to a different representation 
of the quantum state, e.g., taking the Fourier transform switches from the position 
representation to the momentum representation and in this new representation we 
have to consider instead the projection of to momentum space. 

We want to use a method which allows for a bit flexibility, so that we can 
accommodate all possible different representations. To this end we look at how we 
can in general approximate the full Hamiltonian H hy a simpler Hamiltonian Hi in 
a way that the propagation of Lagrangian states initially associated with A can be 
approximated using Hi. The crucial observation which will guide our conditions on 
Hi is that if B{p,q) and A{x) are sufficiently smooth functions and B is the Weyl 
quantisation of B, then 

B(^e*^) = {BaA + 0{h))eT^'^ (71) 

where Bf^[x) = B{S/S{x),x), this is a classical result, see, e.g, [30]. This means that 
the Lagrangian state is concentrated in phase space on A, in particular if B\ — 0, 
then Btp = 0{h). More generally, if B vanishes of order N on A, i.e., {d°'B)\\ — for 
\a\ < N, then 

Bi; = 0{h^) . (72) 

Since a propagated Lagrangian state is concentrated near the propagated 
Lagrangian manifold At — <i>*(A), where $* is the Hamiltonian flow generated by H, we 
expect that if Hi is close to H near At then it will provide an accurate approximation 
for the purpose of propagating states concentrated near At. The conditions for a 
function Hi to be a general flrst order approximation of H near At are 

(i/-i7i(t))|A, =0 , and (Vi/-Vi/i(t))|A, =0 . (73) 

These conditions ensure that the Hamiltonian vector fields generated by H and Hi 
agree on At , hence if we denote by $* the time t map generated by Hi , then 

'i>*|A = $*ilA, (74) 

i.e., the classical dynamics agree when restricted to A. 

Let us now turn to the quantisations of H and Hi. Let U{t) and Ui{t) be the 



time evolution operators generated by H and Hi , respectively, then we expect by ( 74 ) 
that U{t)ip « Ui{t)tp for a Lagrangian state associated with A. To verify this we make 
an Ansatz 

U{t) = Ui{t)V{t) , (75) 
and inserting this into the Schrodinger equation for U{t) gives for V(t) the equation 

ihdtVit) = 5H{t)V{t) , with y(0) ^ I , (76) 



where 



6H{t) = U^{t)[H - Hi{t)]Ui{t) 



(77) 
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By Egorov's theorem, see Appendix A. 3 the operator 6H{t) is in leading order in h 



the quantisation of H — Hi transported along the map $* , i.e., 



SH{t) = [H- Hi{t)] o $* + Ot(fi2 



This implies by |73| that 



5H(t)\A„ = Ot{h^) , and W6H{t)\A„ = Ot{h^) 



(78) 



(79) 



and therefore by ( 72 ) we have for an initial state of the form ip = Aoe r ^° , where Aq 



has bounded derivatives, that 

Wmmw = Ot{h^) . (80) 

Thus the generator of V{t) is small when applied to i/j and therefore Duhamel's 

(81) 



principle, (201, gives 



v{t)^p = ip + Ot{h) , 

hence Ui{t)tp is a good approximation for U{t)ip if is a Lagrangian state. 

This is a generalisation of the standard time dependent WKB method presented 
in Section [2] The choice of Hi corresponding to the results in Section [2] is 



Hi{t; ^, x) H{WS{t, x),x) + W^H{WS{t, x),x)-^ 



(82) 



which we obtain from H{p, q) — H{WS{t, x) 
shifted momentum ^ = p 



one finds using Appendix A. 2 
and so by comparing this wit 



^, x) as a first order approximation in the 

Viq) 



VS'(i, x). For a Hamiltonian of the form H = ^p^ 



that Hi = H{VS{t, x), x) - ihVS{t, x) • V - 'fAS{t, x), 



14| and ( 11 1 we see that in this case we have 



Ui = e 



T{t)e' 



and V{t) = ei^°D{t)e 



(83) 



So how does a caustic affect this relation? In defining Hi we used coordinates (^, x) 
near At defined hy p — VS'(i, x) + £, and q — x, so that At = = 0}, but at a caustic 
this coordinate system becomes singular. In order to repair this we can just use a 
different canonical coordinate system near A(, one which does not become singular, 
and choose for Hi the first order Taylor approximation in the transversal coordinates. 

The only problem we might encounter with this more general version of time 
dependent WKB theory is that the control of the time dependence of the remainder 
term (80) becomes more involved. In Section [2] we were able to take advantage of the 
very explicit version of Egorov's theorem from Appendix A. 3. 2 to reduce this problem 



to the non-contractiveness of the transport map 4iK{t)- For different choices of Hi 
this problem can become much harder and we reserve a more thorough investigation 
for the future. But this question concerns our ability to get explicit error estimates, 
it does not prevent us to use, after suitable testing, the method for explicit numerical 
propagation in concrete problems. 

Now if the initial state V'o is localised in phase space around a point (p, q) G Aq, so 
that the estimate ( 80 ) no longer holds, then it is natural to approximate V{t) by taking 
only the behaviour of its generator SH{t,^,x) near {p,q) into account. By (73 1 we 
have 5H(t,p, q) — and \75H{t,p, q) — 0, hence the simplest nontrivial approximation 
is quadratic, with z = x) and zq = [p, q) we have 



5H{t,x,0 



5H2[t,z)^]^{z 



zo)5H"{t,zo){z-zo) 



(84) 
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The time evolution Mp_q{t) generated by 6H2{t) via ihdtMp,q{t) — SH2{t)Mp,q{t) and 
Mp q{0) = / is a mctaplectic operator since the generator 6H{t) is the quantisation of 
a quadratic function on phase space. Therefore for ip strongly localised around {p, q) 
we expect V{t)il) « Mp^qtp and hence we have arrived now at the more general version 



of (34), which reads 

Uit^j = Ui{t)Mp,q{t)^ + Ot{Vh) (85) 
if il^ is concentrated in phase space at [p, q) as in ( p5|) 



The relation of Mp q(t) to Mq{t) follows from (83), ii Hi is of the form (82) we 
find that 

Mp.qit) = e^^o' Mqit)e-'^'^o'' (86) 

where p = VSo{q) and S'q^-' (x) = p{x — g) + ^ (x — q)SQ{q) {x — q) is the quadratic part 
of the Taylor expansion of Sq around x — q. 

Since Mp^q{t) has a quadratic generator it is the quantisation of a linear symplectic 
map Pp,q{t) on Tp^q{W^ x R") and by construction this map can be expressed in terms 
of the linearisations of the maps $* and (^\. Since V{t) = Ui^{t)U{t) we can view 
V{t) as a quantisation of the map {^\)~^ o $*, and Mp q{t) is then the quantisation 
of the linearisation of that map around p, q, i.e., 

Pp^qit) = (d$*)-M$* . (87) 

Since and $* are identical on A we have that Pp,q{t)\T Aq = I, hence Ppq is a 
shear relative to the tangent space Tp q Aq of the initial Lagrangian manifold at (p, q). 



In case that Hi is given by ( 82 ) the map Pp,q {t) can be described in some more 
detail. Note that since Hi in ( 82 ) contains only a linear term in the momentum we 
find that Hamilton's equation for Hi give 

X = VpH(yS{t,x),x) , (88) 

that means the trajectory of x under $* does not depend on the initial momentum 
^. Hence ^{ maps vertical subspaces into vertical subspaces, i.e., d^\{Vza) C V^jt), 
where for z — {p, q) we set Vz — {{£,, q) G R"}. Therefore the map Pp,q{t) satisfies 

PpAtr'K.Ao ^ I and Pp,,(i)-^(%,)) = (d$*)-i(14*(p,,)) (89) 

and as we show in [Appendix B| this implies that the knowledge of the Lagrangian 
subspace (d<i>*)~^(V$t(p^^)) determines the map Pp^q{t)~^ and hence Pp,q{t) uniquely. 

Let us look at two examples where we can determine the long time limit of 
(d$*)~^(Vit(p^g)) from dynamical conditions. 

(a) If we have a one-dimensional system and d^H/dp^ > 0, then the velocity increases 
with p and hence (d$*)~^(V$t(p,g)) —?■ {{0,q),q G R} for ^ ^ oo. Therefore if 
T(pq)Ao ^ {{0,q),q S R} we can use the result from Appendix B and see that 
for large t the map Pp^q{t) tends to a limit P^q and therefore the corresponding 
metaplectic correction Mp^q(t) tends to a limit, too. 

(b) If the trajectory through {p,q) is hyperboHc and the vertical subspaces Vz 
are transversal to the unstable subspaces along the trajectory, then 

(d$*)-i(y<i, 

'(p,q)) tends to the stable subspace Vp q for large t, and at an 
exponential rate if the hyperbolicity is uniform. So in this situation we find that 
for large t the map Pp^q{t) tends to a limit which is uniquely defined by the 
conditions that J^^lrp.gAo = I and P^q{Vp_q) = Vp^q. Therefore in this situation 
the metaplectic correction tends for large t to a limit, too, and expontially fast if 
the hyperbolicity is uniform. 
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6. Examples 

In order to illustrate how the theory we developed in the previous sections works, 
it is very instructive to look at a couple of simple examples. They will allow us to 
understand in more detail the interplay among the classical dynamics, the choice of 
an initial Lagrangian manifold and the analytical constructions we developed. 



6.1. The free particle 

The first example we look at is the free particle in 1 dimension, with Hamilton operator 
— ^A. We want to propagate a state that is initially localised at {p,q) G R^, where 
p > 0. As initial phase function we choose 

Soix)^pix~q) + ^{x-qf (90) 

where a G K, is a real parameter. The corresponding Lagrangian manifold is 

Ao ^ {{p + ax,q + x) , X eR} (91) 

which is a line through {p,q) with slope a. The corresponding classical Hamilton 
function is H{£^) — and the Haniiltonian flow it generates is given by $*(^,a;) = 
{^,x+tS,). Applying the flow to Aq gives A^ = <&*(Ao) — {{p+ax,q+x+t{p+ax)) , x G 
E,} and by replacing x by x/{l + at) this can be rewritten as 

At ^ {{p + a(t]x,q(t) + x) , X eM.} , where a(t) = ,q(t) = q + tp, (92) 

1 + at 

which is a line through {p,q{t)) — ^'(p, g) with slope a{t). Notice that if a < 0, i.e., 
if the initial line has negative slope, then a{t) — > oo for t — >■ — 1/a, this means that 
the line At turns vertical, hence we have a caustic. But if a > 0, then a(t) < a for all 
t > 0, and no caustics occur. The phase function which generates At and satisfles the 
Hamilton-Jacobi equation with initial condition Sq is given by 

S{t,x)^^-ph + p{x-q{t)) + ^{x~q{t)r . (93) 

To find the transport operator we have to find the map 0a x) — 7rAt$*(7r^J(a;)). 
To this end we notice that tt'^^{x) ~ (VS'o(a;), x) = (p + a(x — q),x), hence 
$*(7r^J(a;)) = {p + a{x — q) , x + t{p + a{x ~ q)) and the final projection just takes the 
second component, therefore 

(l)A{t,x) = {1 + at)x + t{p- aq) . (94) 

Since (j)'f^(t,x) = (1 + at) the condition a > guarantees that (f)\ is non-contracting 
and we get (/)^^(t, a;) — y^^{x — t{p — aq)). Therefore the action of the transport 
operator reads 

Tit)Aix)^^^^A^^^ix^tip-aq))) , (95) 

which for a > is well defined for all t > 0. In order to compute the metaplectic 
correction, which is generated by (49), we notice that A{t,x) = (1 + at)^"^ is 
independent of x, hence we have 

= AW - . (96) 



How do wave packets spread? Time evolution on Ehrenfest time scales. 



17 



Then the metaplectic correction (551 is 



Mg{t)=e-^^^, (97) 

where the time dependent factor in the exponent comes from the integration of the 
time dependence in A{t), (xzp^pp- dt' = j^^- We see that in particular, as we 
observed at the end of Section [5] that the operator tends to a limit for large t, 
Mq{t) Mg{oo) =eT^^. 

We computed all the elements in our extended time dependent WKB propagation 
scheme for the free particle. Let us apply this to an initial state ipoi^) — 
ylo(a;)es[p(^-«)+f with Aq = L^oq for some a e S(R). Then Mg{t)Ao = L^at 

with at = Dq{t)ao and since Dq(t) tends to a limit for large t, we have that at G ^(IR) 
with bounds uniform in i S E,"*". Hence we find 

^P{t,x) = A(i,a;)e*[5P'*+p(^-9W)+4^(:^-9(t))'] ^ (gg) 

with 



Ait, X) = TmmA,[x) = jrrj^^^^^at {j^,J^^^{^ " t(p - .,))) 



(99) 



Since the Hamiltonian is quadratic this expression is actually exact, the remainder 
terms are zero. Notice that the behaviour of the amplitude A{t, x) depends on the 
size of h}^'^{l + at), which determines the Ehrenfest time scale ([s]). li h}/'^{l + at) — > 
the amplitude becomes concentrated, but if h}/'^{l+at) ^ 1 then the amplitude A{t, x) 
is a smooth function and ipi'^) becomes actually a Lagrangian state. 

Let us note that for the special case that our initial state is a Gaussian coherent 
state, 

^o(^) - ^^eslP^'-^^+s^"-")'! , (100) 
(nh) 4 

our result is identical to the one obtained by the standard propagation of coherent 
states, since the Hamiltonian is quadratic. The only difference is that we derived it 
in a different way. But it was still instructive to go through the derivation since it 
highlights a few important points: We can rewrite the state using the phase function 



^0 from (|90| as 

V'(x) = [L5ao](a;)e*^«(^) , with ao(a;) = e'^ , (101) 



and then applying (|98j) and ( |99[ ) gives us the propagated state. But since the initial 
state is independent of a, the final result is independent of the choice of a, too, i.e., of 
the initial Lagrangian sub-manifold used to propagate the state. But the intermediate 
steps depend on a, first of all, we needed to choose a > in order to avoid caustics, 
then the choice of a determined how we split the evolution into transport, implemented 
by the operator T{t), and dispersion, represented by D{t) — Mo{t). The simplest 
choice would have been a = 0, then Ao is horizontal, and actually invariant under the 
classical flow. In this case T{t) is defined by transport along one central trajectory, 
and all the dispersion is described by D{t). If a > 0, then Aq is transversal to the flow 
and the action of T{t) takes actually a family of trajectories into account. Since these 
trajectories move with different speed, this already accounts for part of the dispersion 
of the wavepacket, and the operator M{t) has only to add the dispersion transversal 
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to At- This is reflected by the fact that the coefiicients of A{t) tend to for large t, 
and fast enough so that the hmit 

lim Mo(t) = e'^"^ (102) 

actuafly exist if a > 0. This means that for large t we can replace Afo(t) by its limit 
and still get a good approximation of the propagated state. This is why a > is 
preferable to a = 0, in particular if we extend these construction to more general 
systems. 

In addition, the fact that for a > the coefficients of A{t) tend to for large t 
reflects the property that the map 0a is expanding, this has as well consequences for 
the size of the remainder term if we add higher order terms to the Hamiltonian. To 
study this effect will be the main focus of the next example. 



6.2. Integrable systems 

Let us look now at a more general Hamiltonian which is a function of the momentum 
only, but not necessarily quadratic anymore, so H = H{S^), and we will assume that 
for ^ > we have > and H"{£^) > 0. This is the form a Hamiltonian of an 

integrable system takes in action-angle coordinates. 

The classical dynamics is given by $*(^, x) = {£,,x + tH'{£^)) and by the condition 
H"{C} > the velocity H'{£^) can not decrease with increasing ^. This implies that 
if we take again an initial phase function of the form So{x) = px + ax'^/2 then the 
corresponding Lagrangian submanifold does not develop caustics for i > if a > 0, 
We flnd At = $*(Ao) = {{p + ax, x + tH'{p + ax)) , a; e E,} and therefore 

(l>f^{t,x) = X + tH' [p + ax) , and (f)f^{t,x)~\ + atH"{p + ax) (103) 

so if a > and H" > 0, then (j)\{t,x) is actually expanding. Since the flow is 
now no longer linear. At is no longer a straight line and we do not attempt to find 
an explicit expression for S{t,x), but rather take it as defined by the Hamilton- 
Jacobi equation. But we observe that since At = {(VS'(t, (f>A{t, x)), <j>A{t, x))} we have 
VS'(t, (/)A(i, a;)) = p + ax — V5'o(a;) which is of course a consequence of momentum 
conservation. 

Having the phase function and the map (f>Ait), and therefore the transport 
operator T(t), we have the ingredients of the standard WKB propagation, and 
now we want to compute the metaplectic correction for an initial state localised at 
{p, 0) . Here we have to use the generalised formulation from Section [5] since the 
Hamiltonian is no longer a sum of a quadratic kinetic energy term and a potential. The 
transport operator is the quantisation of the map $*(^,2;) = {[(j)'j^{t, x)^]~^£^, (l)\{t, x)) 
and so the operator SH has symbol SH{t,^,x) — (5ij(°)(t, $* (^, a;)) + Ot{h^) where 
5H(°^ (t, ^, x) = H{VS{t, x) + - H{VS{t, x))~^- VH[VS{t, x)) and the Ot (h^) term 
comes from Egorov's theorem. Therefore using W S{t,<j)\{t,x)) — VS'o(x) we find 

dH{t, x) =H{VoSix) + [Mt, ^)^]-'0 

- H{VSoix)) - [(^'^{t, x)^]-'^ ■ VH{VSoix)) + Ot{h^) . 
Since the point [p, 0) corresponds to ^ = 0, a; = we find using ( 103[ ) that the generator 
of the metaplectic correction is SH2{t) = — ^ [i+atw'\p))^ ^ ^'^'^ hence 

ifi g"(p)t A 

Mo{t) = e ^ i+°*«"<f) . (105) 
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We observe the same phenomenon as in the free case, for a > the large t behaviour 
has a limit. But now this is as well reflected in the behaviour of the remainder term 
which is determined by the size of 

SH{t, \/^, Vhx) - 6H2{t, Vh^, Vhx) = 0{h^/'^[<f)'^{t,x)<f)'^{t,x)^-^) , (106) 

and since the time integral of this quantity is bounded we get for a coherent state ipo 
concentrated at (p, 0) that 

||?7(i)Vo - f/i(t)M(i)Vo|| < C^/h . (107) 

Here we made the assumption that we can control as well the remainder term Otih?) 

from the application of Egorov's theorem, which seems likely since the transport 
operator has still a simple form. But we leave a detailed investigation to a future 
publication. 

This is in contrast to the standard coherent state propagation which can only 
work if \t\\/h <C 1, since Te = ^/Vh for an integrable system. The physical reason 
for this is that for t^/h ~ 1 the state Ui{t)M{t)tpo is no longer a coherent state but 
has become a Lagrangian state, i.e., there is a qualitative change in the nature of the 
state at the Ehrenfest time. 



6.3. A parabolic harrier 

Lot us now look at a case with a hyperbolic trajectory, lot H{S^,x) = — ^x"^, this 
Hamiltonian describes the motion in a parabolic barrier if Vb > 0. The classical flow 
is given by 

^\^,x) = (cosh(At)C + Asinh(At)a;,A"^sinh(At)^ + cosh(Ai)a;) (108) 

where A := ^/Vo > is the Liapunov exponent. The origin is a hyperbolic flxed point 
with stable and unstable manifolds given by 

V = {{-Xx,x),xG'R} , V = {{Xx,x) ,xGR"'} . (109) 

Let us first choose an initial state localised on the fixed point, i.e., at (0,0). An 
initial phase function So{x) = a^x^/2 corresponds to Aq = {(ax,x) , x S E,} and the 
transported manifold At = $*(Ao) can be written as 

r/ / N N . , /X cosh(Af)a + Asinh(Af) 

At = {(a(t)x,x),xeR} , with a(t) = , , .\ , ttvt • (HO) 

■■^ ^ ^ ^ -* A-i sinh(Ai)a + cosh(At) ^ ^ 

The corresponding phase function is S{t,x) = a{t)x'^/2 and we flnd furthermore 

0A(t, x) = (aA^^ sinh(Ai) + cosh(Ai))a; , (111) 

and therefore 

Ait) = (aX-^ sinh(At) + cosh(Af))-2A . (112) 

The manifold At does not develop a caustics if a > —A. The case a = — A is 
special since then = V, i.e., A^ = Aq is invariant and equal to the stable manifold. 
In this case the manifold is contracting, but for all a > —A the manifold is non- 
contracting. The case a = A is as well special since then A( = for all t, so the 
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initial manifold is the unstable manifold, in all other cases for large t. These 

different cases arc as well reflected in the behaviour of T{t) and M{t). For a = —A 
we have 4'a{x) = e~^*x and hence 

T{t)A{x)=e^*/^A{e^*x) and M(i)=e^'^^. (113) 

The transport operator T{t) is squeezing at an exponential rate, and the metaplectic 
correction has to make up for this at an exponential rate, too. In contrast, if a > —A 
we have 

qA-1 sinh(Ai) + cosh(At) = ^^e^* + 0(0-^') (114) 

and so T{t)A(x) = [(^^(t)]~^/^A((/)^^(t,a;)) is stretching at an exponential rate. The 
metaplectic correction is given by 

M(i)=e^'(^+°(^"''^1^ (115) 

and tends exponentially fast to a limit. This dichotomy between the cases a = — A 
and a > —A is analogous the the one between a = and a > we found for the 
free particle. For a > — A the metaplectic correction saturates in time, whereas for 
a = —A its contribution grows exponentially. Since the Hamiltonian is quadratic 
both cases give exact solutions, and the different initial manifolds only correspond to 
different ways to split the time evolution. But as in the integrable case, if we move 
to a perturbation the contracting case a = —A gives remainder terms which blow up, 
whereas the expanding cases a > — A give remainder terms which remain bounded by 
Vh independent of time. 

So far we have looked at an initial state which is concentrated on top of 
the barrier, let us now look at a initial state localised at {p, q) wit g < and 
p > 0, hence an incoming state. The question of interest is then if this state gets 
transmitted or reflected, and how one can describe the transition between reflection 
and transmission uniformly if one changes {p,q). We chose an initial phase function 
So(a;) = p{x — q) + a{x — qf jl, which again for a > —A gives an expanding initial 
manifold Aq. For simplicity we will choose a = A, i.e., Aq = (p, <z) + V^", the general 
case a > —A can be treated similarly. From computing $*(V5o(a;), x) we find 

Mt,x)=e^\x-q) + q{t) , (116) 

where q{t) = X~^psm]i{Xt) + qcosh{Xt), and therefore 

T{t)L,a{x) = ^^a^-^[x-q{t)]y M(i)=e5^^^. (117) 

For the phase function we find 

S{t,x) = C{t)+pit){x-qit)) + \{x-q{t)f/2, where {p{t), q{t)) = 'i>\p, q) , (118) 

and C{t) = /q (pg - H{p, q)) dt. We see that M{t) again tends with exponential speed 
to a limit M°° = e^A^. The behaviour of T{t)Lqa depends crucially on e~'^*/\/h, 
for small time this parameter is large, hence the state is localised, but for 



(119) 
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i.e., on Ehrenfest time scales, the function T{t)Lqa is no longer strongly localised and 
the transported wave packet is actually a Lagrangian state associated with 

At = $t(Ao) = (p(i),g(t)) + ^" . (120) 
Where this state moves depends on q{t), and for long times we have 

q{t)^^^e'' + 0{e-'') (121) 

so the behavior oiq{t) depends on the value of p+Ag. The cs&c p + \q = corresponds 
to initial conditions {p,q) in the stable manifold of the fixed point at (0,0), and 
so the trajectory runs into the fixed point and the Lagrangian state is for long times 
located on the unstable manifold V". The case p + \q > gives a wave packet 
which is transmitted over the barrier and the case p + Xq < corresponds to a wave 
packet which is reflected. The strength of our approach is that we can describe the 



transition between these different cases in a uniform way given by (117). We studied 
a simple example here, but the method works for general potential barriers and allows 
to describe the transition between reflection and transmission of time dependent wave 
packets in a uniform way. 

Such processes are of great importance in the theory of chemical reactions and 
we expect that our method could be of great use in the theoretical and quantitative 
description of time resolved chemical reactions which are experimentally accessible in 



femto and atto chemistry 34 35 . A first step will be to analyse the transport of wave 



packets over barriers, or more generally through bottlenecks in phase space, using the 



method of Quantum Normal Forms recently developed in 36 37 . 



7. Transport along a hyperbolic trajectory and its unstable manifold 

In this section we want to discuss the case that the initial coherent state is concentrated 
in a point z = {p,q) on a hyperbolic trajectory z{t) — <i>*(z). As a test system we 
choose the kicked harmonic oscillator 

^ oo 

H{t,p,q,t) = -{p^+q^)+Kcosq ^ S{t - n) , (122) 

n— — oo 

where K is the chaoticity parameter which we will choose to be if = 2. This system 
has as well been used in 24 . In part (a) of Figure[l]we show a portrait of the classical 



phase space, the system has a an unstable fixed point at the origin and we display the 
unstable manifold. The Liapunov exponent of the unstable fixed point is A = 0.83... 
and we use h = 0.0008, so the Ehrenfest time is 

Te = 4.29... (123) 



As initial state we choose a Gaussian wavepacket centred on the fixed point at 
the origin: 

In Fig. [l] (b) we display the Wigner function of the time evolved state at t = 4 which 
can be seen to be stretched along the unstable manifold. As expected at the Ehrenfest 
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q 



Figure 1. In (a) we show a phase space portrait of the KHO dynamics, | |122| l, 
for K = 2. The origin is a hyperbolic fixed point, and we displayed as well 
the corresponding unstable manifold. In (b) we plot the Wigner function of 
the evolved state ( |124| at t = 4^ , which is of the order of the Ehrenfest time 
= 4.29 .... It evolves along the unstable manifold and starts resembling a 
WKB state. For comparison in the upper left corner the Wignerfunction at t = 
is shown, note that the vertical axis has been rescaled compared to (a). 



time the state becomes a WKB type state associated with the mistable manifold. For 
comparison the real part of the time evolved wave function iplt) is plotted in part (c) 
of Fig. [2] 

The prediction of the theory in Section [3] is that, after a suitable metaplectic 
correction of the initial amplitude, the state can be propagated using standard time 
dependent WKB approximation . This means that, if we decompose the initial state 
as 

V'o(a:)-(ioa)(a;)e*^°(^)/\ (125) 



then the evolved state is by (19) and (34) given by 



^{t, x) = [T{t) D{t) Loao] (a:)e^^(*'-)/'' « [T{t) Mo{t) L„a] (x)e^^(*'-)/^ . (126) 

Here S{t, x) is the solution of the Hamilton- Jacobi equation satisfying S{0, x) — So{x). 
The operators T{t), D{t) and Mo(t) also depend on the choice of S'o('7)- Even if 
according to Section |4.4| this choice is, to some extent, arbitrary, the efhciency of 
the method may be sensitive to Sq{x). We shall only consider the quadratic case 
Sq{x) — ^^f^x^, corresponding to a linear Lagrangian manifold p = tan(0) g. We 
shall see that there is a wide range of initial manifolds that work well. 

Note that Eq. ( |126 ) is equivalent to saying that the sequence of operations 



T*{t) U{t) e-'^(*)/'^l = D{t)Lo ao « M{t)Lo oq 



(127) 



must produce approximately a Gaussian state centred at 0, for a Gaussian initial 
amplitude ao{x). This is the first test we shall carry out, with the choice So — 
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Figure 2. In the top panel we test l |127| l: The initial state | |124[ l is propagated 
forward using the full propagator to t = 4 and propagated backwards using 
just time dependent WKB, the result (symbol) is compared with the metaplectic 
approximation (full line): (a) amplitude, (b) phase derivative. We also show the 
initial state for reference (dotted lines in (a)). In the bottom panel we compare the 
metaplectically extended WKB scheme with exact quantum propagation (circles). 
The right panel (d) is a zoom of the tail of the wavefunction in (c) to show that 
agreement extends down to this scale and into the most nonlinear region 



corresponding to the initial Lagrangian manifold p = 0. First we propagate the state 
( 125 ) exactly during some time, and then propagate it back using time dependent 
WKB. Panel (a) and (b) in Figure [2] shows the comparison of the exact result D{t) LqUq 
with the metaplectic approximation Mf){t) Loa^. In the case of the modulus the 
agreement between the exact DLoUq and MoLgao is perfect within visual resolution. 
There is a small difference between the phases, but this occurs in a region where the 
amplitude is very small. We show as well a direct comparison between the exact and 
the metaplectically extended WKB propagated states in part (c) and (d) of Fig. [2] 
The agreement is excellent. 

We will now investigate how robust the method is with respect to the changing the 
choice of the (linear) initial WKB manifold. We have changed the slope 6* in a range 
of almost 90 degrees, from 6 = — 0.307r/2 to 6* = 0.657r/2. The comparison between 
exact and metaplectically extended WKB at < = 4 are shown in Fig. [3] Except for 
very slight differences (blue triangles) the agreement is still excellent for these "large" 
slopes. 
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Figure 3. We compare the effect of different initial Lagrangian manifolds 
p = tanS g on the propagation, we considered for 8 the following fractions of ■it/2: 
-0.30, 0.35, 0.65. On the left we display \DLoao\{q) (circles) and the metaplectic 
approximations |MoLo''ol(<7) (lines). On the right we show the evolved state at 
t = 4~ (real part). We compare metaplectically extended WKB (symbols: black 
dots and blue triangles) with exact quantum propagation (line), two slopes where 
considered: 0.657r/2 and — 0.307r/2 

As we increase the slope, the amplitude concentrates in narrower regions. This 
happens because, the larger the slope, the map induced in the g-coordinate by the 
WKB manifold is more expansive, and T* is more contracting. 

We tested the theory as well for different values for h and found good agreement 
with the expected behaviour (not shown). 

8. Conclusions 

In this paper we derived an extension of the standard time dependent WKB method 
which can be applied to highly localised states like coherent states. It allows to describe 
in a uniform way the transition in time from a semiclassically highly localised coherent 
state to an delocalised Lagrangian state which takes place at the Ehrenfest time. 

The main idea on which this extension of time dependent WKB theory is built 
is an exact decomposition of the time evolution of an initial state of the form 
i'oix) = AQ{x)e'K^'>'-^\ where Sq is real valued, into several parts 

m = (T(<)i?(t)Ao)e^^W . (128) 

Here S{t) is a solution of the Hamilton Jacobi equation and is hence related to the 
transport of the Lagrangian manifold Aq generated by 5*0 through phase space. The 
unitary operator T{t) transports functions in position space along the projections 
of the phase space trajectories emanating from Ag. Finally D{t) is the propagator 
generated by the time dependent Hamiltonian —H'^T*{t)AT{t)/2. S{t) and T{t) are 
defined purely in terms of transport along classical trajectories and D{t) takes into 
account the dispersive effect of quantum mechanics. 

The standard time dependent WKB approximation is obtained by approximating 
D{t) /, this works fine if the amplitude Aq is sufHcient flat, i.e., has bounded 
derivatives. But for a coherent state Aq is strongly localised around a point x = q, 
then it is more natural to freeze the coefficients of the generator of D{t) at x — q, the 
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resulting operator is a metaplectic operator Mq(t) whose action on functions can be 
computed quite easily. The main result of this paper is thus 

m « {T{t)Mg{t)Ao)e-^-'^'^ . (129) 

which is valid for amplitudes Aq which are strongly localised around x — q. In order 
to justify the validity of this approximation for long times, particularly on Ehrenfest 
time scales, we had to analyse the underlying classical dynamics more carefully. We 
introduced a non-contraction condition on the position space trajectories emanating 
from a neighbourhood of the initial state with momenta p — VS{x)^ and if this 
condition holds our propagation scheme is effective. The non-contraction condition 
excludes caustics, but we have a large freedom in the choice of the initial phase function 
S{x) which allows in many cases to avoid caustics, at least until the state becomes 
delocalised. 

For times shorter than the Ehrenfest time our scheme reproduces the standard 
coherent state propagation results based on a Taylor expansion of the Hamiltonian 
around the centre trajectory. But for times of the order of the Ehrenfest time we 
find that the state becomes extended and a Lagrangian state, which in the chaotic 
case is supported by the unstable manifold of the centre trajectory. From that time 



onwards standard time dependent WKB theory applies, as has been observed in 24 
In particular if the system is hyperbolic and mixing one can apply the results from 28 
to conclude the state becomes equidistributed after the Ehrenfest time. 

In order to extend the results to more general Hamiltonians, and to be able 
to include caustics, we noticed that standard time dependent WKB approximation 
can be viewed as the exact quantum time evolution generated by a quantisation of 
a first order Taylor approximation of the Hamilton function around the Lagrangian 
manifold associated with the time evolved WKB state. Based on this insight we 
could now choose different first order approximations which remained valid for general 
Hamiltonians and at caustics. The price one has to pay is a more complicated and 
less explicit formalism. In addition all the previous results were in principle rigorous, 
although we refrained from stating them in the form of theorems, but here we have 
to rely on an assumption that certain special cases of Egorov's theorem remain valid 
on Ehrenfest time scales. But we get a further benefit from this more general way 
to look at the time dependent WKB approximation, we can show that the classical 
map associated with the metaplectic correction Mqit) is a shear map on phase space 
and that it furthermore in many cases converges to a limit for large t which can be 
determined from simple geometric considerations. This implies that in many cases 

lim M„(t) = M(°°) (130) 
and we can determine Mq°°^ up to a phase from simple geometric considerations. If 



the centre trajectory is hyperbolic then the limit in ( 130 ) is reached exponentially fast. 

We illustrated and tested the theory with several examples. For the free particle, 
and more generally, for integrable systems, the dynamics can be computed quite 
explicitly. The Ehrenfest time is of order Te ^ h'^^^^ and we see a qualitative 
transition in the nature of the state from a localised coherent state to an extended 
WKB state at this time scale. Similarly for a parabolic barrier the dynamics can be 
computed explicitly and the formalism developed in this works gives explicit formulas 
which describe in a uniform way the transition from reflection to transmission of a 
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wavepacket when one varies the energy near the critical energy. We then carried 
out detailed numerical tests on the kicked harmonic oscillator for an initial state 
localised on a hyperbolic fixed point. These showed impressive agreement between the 
metaplectic extension of the WKB method and exact quantum propagation. These 
tests illustrated as well the simplicity of the method; the metaplectic correction is in 
fact a very simple operator, and if one has implemented the time dependent WKB 
method then adding the metaplectic correction is easy and allows at once to propagate 
a much larger class of states. 

There are many areas where the results from this paper should be of interest. 
In many physical system the Ehrenfest time is for a realistic set of parameters 
actually quite short. E.g., in quantum billiards coherent states spread out after 



a few bounces, 22 23 , and the metaplectic extension of WKB should be able to 
describe this transition. The scattering of a wavepacket of a barrier near the critical 
energy is another example, the Ehrenfest time is the time when the wavepacket reaches 
the barrier, i.e., the time where the physically interesting processes start to happen. 
Metaplectically extended WKB allows to describe this process explicitly and uniformly 
in h and t. A large class of chemical reactions is described in the framework of 
Transition State Theory by the crossing of a barrier on a high dimensional energy 
surface and modern experimental methods in atto and femto chemistry allow to study 



the dynamics of such chemical reactions with impressive precision, 34 35 . The 
metaplectic extension to WKB together with the normal form approach to Transition 
state theory, 1 36 , 37^ , should allow to give an efficient theoretical description of chemical 
processes on such time scales. 



Appendix A. Wigner Weyl correspondence 

In this appendix we collect some material on Weyl quantisation which we use in the 
main part. Most of it is a development of standard material and only the result on 



Egorov's theorem seems to be new. References for the material we present are 29 30 



Appendix A.l. Weyl quantisation 

Let W{(,,x) := es(^'J+^P) be the Weyl operators which represent translations in phase 
space, then we can define for a function A{p, q) on phase space an operator, its Weyl 
quantisation^ by 

A = jj TA{tx)W{i,x)^^ (A.l) 

where TA{S^,x) ~ JJ ^(p, g)c^s-[p5+9=^l dpdq denotes the Fourier transform of A. The 
function A{p, q) is called the Weyl symbol of A, and the properties of A can often be 
determined from properties of A, e.g., if A G S"(E" x R"), then A : S{M") S"(E"). 
li A = \tp){'4'\ is the projection onto a state ip, then A{£^,x) is proportional to the 
Wigner function of the state tp. 

The product of operators can be expressed in terms of the symbols, one has 
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AB = ^tJB where 

AiBip, q) ^A{p, g)e*[^''-^«~^«-^''lB(p, q) 

=A{p,q)B{p,q) + ^^{A,B}{p,q) (A.2) 

- \a{p, • - V, • ^pfB{p, q) + --- . 

Here {A,B}{p,q) denote the Poisson bracket, and the arrows over the derivatives 
indicate if they act on the function on the left or on the right. Suitable conditions on 
the functions A and B under which this expansion holds can be found in f29|. 

Appendix A.2. Dynamics generated by first order operators 

Let K{t,x,(,) — X{t,x) ■ where X{t,x) is a time dependent vector field, we are 
interested in the time evolution T{t, s) generated by 

K{t) = -ihX{t, x)-V - ^ VX(t, x) . (A.3) 

Let (l>{t,s,x) be the family of maps generated by the vector field X{t,x)^ i.e., 

dt(l){t,s]x) = X{t,(l){t,s;x)) , (j){t,t;x) ^ I , (A.4) 

then 

{T{t,s)A){x) = [det(j)'{t,s;x)]'^^^A{(j)~\t,s;x)) . (A.5) 
Appendix A.3. Egorov's theorem 

Egorov's theorem is one way to formulate the correspondence principle, it gives a 
general relation between classical and quantum dynamics. Let Hit) be a real valued 
smooth phase space function, H{t) its Weyl quantisation and $* and U{t) be the 
classical and quantum time evolution generated by H{t) and H{t)^ respectively. Then 
Egorov's theorem states that 

U*{t)AU{t)^At + Ot{h^), (A.6) 

where At = A o $* , and A and H have to satisfy some conditions on their smoothness 



and growth at infinity, see 27 . Since At is the classical time evolution of A this means 
that quantum and classical evolution are close for small h. The remainder term Otih^) 
does depend on time, and the best general estimates are of the form Ot{h^) ^ S^e'"* 
for some constant F > which depends on H. We want to discuss two cases in which 
one has better control over the remainder. 

One approach to Egorov's theorem is based on writing Heisenberg's equation of 
motion for A{t) — U* {t)AU{t), i.e., ihdtA{t) = [H,A], in terms of the symbols using 



(A.2 ) which gives 

dtA{p, q) = Ia{p, q) sm{h[^p -^q-^q- ^p]/2)H{p, q) 

= {A, H){p, q) + —A{p, • - V, • ^pfH{p, q) + ■ ■ ■ . 

The leading order term is just the Liouville equation and gives A k Ao ^ the main 
problem is then to control the higher order terms. 
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Appendix A. 3.1. Quadratic Hamiltonians and Metaplectic operators If H{p,q) is a 
quadratic function of p and q, with possibly time dependent coefficients, tlien the h 



expansion in (A. 7) terminates after the leading term, and hence the evolution equation 



for A is just the classical Liouville equation. So in this case 

U*{t)Au{t)=At , (A.8) 

where At = Ao^*, and one says Egorov's theorem is exact. The corresponding classical 
maps are linear, hence the operators U{t) for all quadratic H form a quantisation 
of the symplectic group, which turns out to be a double cover of the symplectic group 
called the metaplectic group. The operators U{t) are often referred to as metaplectic 
operators [3 31 32 . 



Appendix A. 3. 2. Conjugation by a flow The second case we need is that U(t) = 
T{t,0), i.e., H{t) is linear in p and given by H{t,p,q) ~ K{t,p,q) = X{t,q) - p. The 
classical map ^\{p, q) generated by H is given by the solutions to 

iit) - ~VgK{t, , Ht) = VpX(t, m,^it)) (A.9) 

with initial conditions ^(0) — p and a;(0) — q. We can express ^{ in terms of the map 
(j){t,0,x) as 

^) = ([0'(O, t, x)]te, 0(t, 0, x)) (A.IO) 

where [4>'{0, t, x)]^ is the transpose of the inverse of the matrix (t)'{t, 0, x). 

Let us consider first the case that A is linear in p, i.e., A = b(q) -p for some vector 



valued function b{q), then only the leading order term in (A.7) is non-zero, and hence 



T*{t,Q)AT{t,0) = At, with A{t,q,p)=b{cl>{t,0;q))-[(l>'{Q,t,x)]^p , (A.ll) 

without any remainder terms, hence Egorov is exact again. Using this result we can 
discuss the case we will need, namely the case that A is a second order differential 
operator of the form 

i = (6(g)V) • S(g)V (A.12) 
So that T*At = T*B{q)VT ■ T*B{q)VT and if we denote the rows of B by hj{q) we 
can use the previous result ( |A.lip and ( |A.2[ ) to obtain T*(i, 0)Ar(t, 0) = A{t) with 



A{t) = HHt,0;q)) ■ W{0,t,x)]^p)i{b,{cj,{t,0-q)) ■ [<^'(0,i,x)]V) 

i 

= p ■ cj^'iO, t, x){B^B)m, 0; qm'{Q, t, x)\^p 

+ y E (^^ • ['^'(0' ^' (^P • ^,){^, ■ ^p) {M<l>{t, 0; q)) ■ [cl^'iO, t, x)]^p) 

(A.13) 



where we have used (A. 2 ) . This can be further simplified, but we will restrict ourselves 



to the case B — I, hence bi — and then wc find 

^(t,p,g)=p-0'(O,i,x)[0'(O,t,a;)]tp+^^Tr[0:'(O,i,q)]2 , (A.14) 

which is the symbol of T* (i, 0)Ar(i, 0). The term of order contains second 
derivatives of the inverse of 0; q) and so if 0; q) is non-contracting, then these 
derivatives stay bounded. 
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Appendix B. A Lemma on Symplectic Maps 

The results in this Appendix are used at the end of Section [5] to show that the 
classical map associated with the metaplectic correction is uniquely determined by 
the dynamics of the tangent space to the initial Lagrangian submanifold and the 
vertical subspace. 

Lemma Appendix B.l. Let Li,L2 C R" x R" be Lagrangian subspaces with 
Li n L2 — {0} and let L be another Lagrangian subspace with L H Li — {0}. Then 
there exist a unique linear symplectic map T with Tj^^ = / and T{L2) = L. 

To show this we will use mainly two standard facts from linear symplectic 
geometry, see e.g. 



30 



(a) By Darboux' Theorem there exist symplectic coordinates {v^w) E R" x R" such 
that Li = {w — 0} and ^2 = {"^ = 0}. 

(b) If L3 is a Lagrangian subspace with L3 n Li = {0} then there exists a unique 
symmetric matrix A such that L3 — {v = Aw}. 

Let us now prove the lemma. The condition T|l^ = / implies that in the 



coordinates from (a) the matrix representing T is of the form Mt = j^j where 
A, B are nxn matrices. Now this matrix must be as well symplectic, i.e., M^UMt — ^ 
^ and this gives B = L and A = A*'. Therefore 

r(L2) = {v = Aw} , (B.l) 



but by (6) the condition T{L2) = L determines then A uniquely. So Mt 
is the unique shear relative to Li which maps L2 to L. 
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